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1 



In this paper we consider the Hamiltonian formulation of the equa- 
tions of incompressible ideal fluid flow from the point of view of optimal 
control theory. The equations are compared to the finite symmetric 
rigid body equations analyzed earlier by the authors. We discuss var- 
ious aspects of the Hamiltonian structure of the Euler equations and 
show in particular that the optimal control approach leads to a stan- 
dard formulation of the Euler equations - the so-called impulse equa- 
tions in their Lagrangian form. We discuss various other aspects of the 
Euler equations from a pedagogical point of view. We show that the 
Hamiltonian in the maximum principle is given by the pairing of the 
Eulerian impulse density with the velocity. We provide a comparative 
discussion of the flow equations in their Eulerian and Lagrangian form 
and describe how these forms occur naturally in the context of optimal 
control. We demonstrate that the extremal equations corresponding 
to the optimal control problem for the flow have a natural canonical 
symplectic structure. 

1 Introduction 

In this paper we consider the Hamiltonian formulation of the equations of 
incompressible ideal fluid flow from the point of view of optimal control 
theory. Our goal in this work is to compare the fluid equation arising in 
this fashion with the symmetric generalized rigid body equations derived 
in Bloch and Crouch [1996], Bloch, Brockett and Crouch [1997] and Bloch, 
Crouch, Marsden and Ratiu [1998]. In these papers we showed that a natural 
control approach led to a form of the rigid body equations on SO(n) x 
SO(n) rather than on T* SO(n). In contrast here we show that the optimal 
control approach leads to a standard formulation of the Euler equations - 
the so-called impulse equations in their Lagrangian form. A nice survey 
of the impulse equations in the various forms can be found in Russo and 
Smereka [1999] (see also Kuz'min [1983], Osledets [1989] and Haddocks and 
Pego [1995] for related work). Lie-Poisson reduction is an important tool 
in rationalizing many of these approaches (see e.g. Marsden and Weinstein 
[1984]). In particular, the Hodge projection to the divergence free vector 
fields, is a Poisson map since it may be regarded as the dual of the natural 
inclusion (this is a standard result; see Marsden and Ratiu [1994]). Thus, 
the Hodge projection naturally takes the unconstrained Poisson system to 
the constrained one. Russo and Smereka and others concern themselves with 
the numerical aspects of these equations, something we do not consider here. 

These impulse equations are not symmetric in the same sense as the 
symmetric rigid body equations - i.e. we do not obtain two symmetric 



equations evolving on two copies of the difFeomorphism group, but it is 

possible to get a more symmetric formulation which we intend to discuss in 
a forthcoming publication and that we mention briefly in the conclusions. 

In the remainder of the introduction we recall the standard and sym- 
metric rigid body equations. 

We recall from Manakov [1976] and Ratiu [1980] that the left invariant 
generalized rigid body equations on SO(n) may be written as 

Q = Qn 

M = [M, n] , (RBn) 

where Q G SO(n) denotes the configuration space variable (the attitude of 
the body), $7 = Q~^Q G so{n) is the body angular velocity, and 

M := J(0) = m + VLke so{n) 

is the body angular momentum. Here J : so{n) — > so{n) is the symmetric 
(with respect to the above inner product) positive definite operator defined 

by 

j(o) = m + fiA, 

where A is a diagonal matrix satisfying Aj + Aj > for all i ^ j. For n = 3 
the elements of Aj are related to the standard diagonal moment of inertia 
tensor / by /i = A2 + A3, /2 = A3 + Ai, Is = Ai + A2. 

The equations M = [M, fi] are readily checked to be the Euler-Poincare 
equations on so{n) for the Lagrangian 

The left invariant symmetric rigid body system is given by the first order 
equations 

Q = Qn 
p = pn 

where O is regarded as a function of Q and P via the equations 

n := J-^{M) G so{n) and M := P - P^Q. 

These equations can be derived from the following optimal control prob- 
lem: 



(1.1) 



Definition 1 Let T > 0, Qo,Qt ^ SO(n) be given and fixed. Let the rigid 
body optimal control problem be given by 

min i f{U,J{U))dt (1.2) 

subject to the constraint on U that there be a curve Q{t) G SO(n) such that 
Q = QU Q(0) = go, Q{T) = Qt. (1.3) 



Proposition 2 The rigid body optimal control problem || has extremal evo- 
lution equations where P is the costate vector given by the maximum 
principle. 

The optimal controls in this case are given by 

U = J-\Q'^P-P^Q). (1.4) 

In §|3| we derive the impulse equations for fluid flow from the optimal 
control point of view. 



2 Inviscid, Incompressible, Fluid Flow 

In this section we introduce the usual dynamics for inviscid, incompressible 
fluid flow, impulse density and the vorticity dynamics. The basic equations 
we consider are: 

dv 

— — \- {v ■ gvad)v = — gradp; divv = (2-1) 



x £ Q; v = v{x,t), p = p{x,t). 

We assume, for simplicty only that the flow is in all of space or in a periodic 
box so we do not need to deal with boundary conditions. This is not an 
essential restriction. 

Here, v is the fluid velocity and p is the pressure. We introduce the 
impulse density z, 

z = V + grad/c. (2.2) 



where k is an arbitrary scalar fleld, k = k{x,t). Notice that the preceding 
equation gives the (Ilelmholtz)-Hodge decomposition of z. In other words. 



the projection of z to v is the Hodge projection of z. We return to this 
important remark in the conclusions to gain deeper insight into what is 
going on with the calculations to follow. 
Take the time derivative of ( |2.2| ) to get 



where 



dz 

— — f X curlz = gradA, divv = (2-3) 



A dk 1 

A is called the gauge. Any choice of gauge is possible, but to be concrete, 
we consider the "geometric gauge" A = —v ■ z. 
With this choice 

— + (t" • grad)z + (gradv)^ z = 0, divv = 

and k is now fixed by the equation 

dk 1 

— = p V ■ V. 

2 

Let z = '^^Zidxi{= z ■ dx) be the one form corresponding to z. Then one 
can easily show that 

£vZ = + z) ■ dx. 
Thus the impulse density is governed by: 

^ + £vz = 0, div V = 0. (2.4) 

Hence by applying the exterior differential operator we obtain 

d 

Trdz + C^dz = 0, div v = 0. 
ot 



Lemma 3 w = curl z = curl v satisfies the vorticity equation: 



(2.5) 



Proof dz = Y^ - dzi A dxi = Y^j^j{curlz)ijdxi (8) dxj, where ab = a x b. We 
may now compute: 

A 



Cydz = ^v, curlz 



dxi dxj . 



We now quickly review the two coordinate systems associated with the 
fluid system. We denote the Lagrange or material variables by Xi and the 
Euler or spatial variables by Xi, and set 

Xi = (f)i{X,t), 1 <i <3. 

We assume ^ : J7 ^ is a volume preserving difFeomorphism, with Jacobian 

equal to unity, = 1. 

Let v{x,t) = spatial velocity, so that 

^=Vi{x,t). 

Thus V{X, t) = v{(p{X, t)) is the material velocity. Hence 

^{X,t)=v{cl>{X,t)) 

or 

^ = vocl>. (2.6) 

We note the "right invariance" of this system and its evolution on the 
"Group" G = DifF^Qj(O) of volume preserving diffeomorphisms of Q,. Setting 

{a,b) = / (x,t)b(x,t)dx, 

and using = 1, we obtain the identity 

(a o ^, 6) = (a, b o (f)~^). 
With this introduction we may introduce the total vorticity equations: 
d(f) dw r 1 1- 

'dt^'"° = [w,v\:dwv = 0. (2.7) 

It is interesting to compare these equations with the right invariant Euler 
equations for the rigid body: 

Q = 1)Q; M= [O, M] (2.8) 



[n, M] = VLM - MQ. 



[M, $7] interpreted 
as vector fields 



These Euler equations are Hamiltonian on T*S0(3), with the canonical sym- 
plectic structure. The equivalent statements about (^) have been well 
studied, (see references in the introduction). However, the derivation of the 



symmetric version as in (1.1) provides our motivation for this new study. 



3 Optimal Control Problem 

In this section we introduce an optimal control problem and discuss the 
corresponding extremals. The problem can be posed as: 

1 r 
mm- / {v,v)dt 

-^i-) ^ Jo 

subject to: 

divf = 0; ~ ^ ° ^ i"^-^) 

and 

</)(X, 0) = MX), HX, T) = MX) fixed, 

and, for flow in all of space, suitable conditions at infinity. 

This optimal control problem is of course identical to the standard Hamil- 
ton principle for ideal fluid mechanics. However our goal here is to analyze 
it from the point of view of the Pontryagin maximum principle. 

We solve this problem by introducing Lagrange multipliers and the cost 

T// 1 



J{v,(l),Tr,k) = I [ {Tr,vo(l)- — \ - -{v,v) + {k,dwv) ]dt 



The problem (3^) may be recast as: min J, subject to divv = 0, = f o(/), 
and boundary conditions. 

We may prove the following result: 



Theorem 4 The extremals of problem (|3.l| ) are given by 
dn rp d(l) 

— = - {V^0(f)) TT, — =VO(f> (3.2) 

V =TT o (f>^^ — grad k, div v = 0. ■ 



Sketch Proof 



d 

6J = {{'K.dv o (j) + {v^ o (j))5(t) — —5<j)) — {v,5v) 
Jo 

div dv))dt 

T 

((vr, 6v o (p) — {v, 6v) + {k, div 5v))dt 
d 

(vr, {v^ o (/))(5(/) - —5(t))dt 
at 

Noting that 5v{oo,t) = 54>{x,ti) = 5(j){x,T) = 0, we obtain 

6 J = / {tt o (p ^ — V — grad k, 5v)dt 
Jo 

+ / {(v^o(l))'^'K + ^,5<t))dt. 



The system (^]^) follows immediately. | 



We note that the system ( |3.2D should be interpreted in terms of Lagrange 
and Euler variables in the form 

^(X,t) =- (^^{ct>{X,t),t)y 7r{X,t), 

^{X,t) =viHX,t)), 

and 

v{x, t) = Ti o (f)~^{x, t) — grad k{x, t). 

We now study the Hamiltonian for the extremal flow. Employing the 
maximum principle we know that the Hamiltonian corresponding to the 
problem (3.1) is 



i^- (vr, (/)) = (vr, w o 0) - v). 

We introduce the vector potential for v 

V = cnvlil); div'(/' = 0, 



(and ip ^ at infinity in all of space) . Thus 

uj = cur\v = curl curl-i/' = —Aip + grad divip = —Aip, 

so 

oj = —Alp; A = Laplacian; . 

Thus tp = Auj where A is an integral operator. From these identities we may 
write: 

H{tt,4>) = {TTO(j)~^^ cur 1-0) — -(f, curl^/^) 
= (curl vr o ■(/;) — -(curl t), -0) 



But 



so 



and 



V = IT o (f> ^ — grad k 



curl V = curl vr • 



H{Tr,(t>) = ^ (curl vro 0-1, V') (3.3) 

= i(7ro0-i,r;) (3.4) 

= \{u^,Auj) (3.5) 

= -(curlvr o 0~i, ^ curlvr o 0~i). (3-6) 

We now compute along extremals (|3 
d 



dt 



Y,Mx,t)dMx,t) = o 



or 

d_ 

dt 



Hence 



Thus the "canonical two form" 



i ijk 



is constant along extremals. We now have the following critical result. 
Lemma 5 Let z = n o 0^^. Along extremals ( |3.2D 

dZ rj. 

+ z^v + z = 0. m 

Thus z = z ■ dx satisfies 

-+£.z = 0, ^=^o0 (3.7) 

It follows that we have recovered the evolution of the impulse density 
of the fluid flow, equation (|2.4| ). Note that H = ^(z, f), where —z ■ v = K 
is the geometric gauge. The following result relates z to the canonical two 
form. 

Lemma 6 4)^"^* ^ ■ dvTj A d(()i = Yli dzi A dxi = dz. ^ 

Substituting the relation 

z = V + grad k 



into the system (3^) recovers the system where 



^=p-2V.v. (3.8) 

However k is not arbitrary, since it is determined by the extremal system 
(|3ll) . In fact 

div z = div vr o (p^^ = div grad k = Ak. 

Thus, 

A: =^div z = Adw vr o (3.9) 
V =TT o (f)~^ — grad Adiv vr o (p^^ . 

Hence the pressure p is also determined by the flow. 



Lemma 7 By augmenting the cost functional in the optimal control problem 
( |3.lD by a potential function rj 



the extremal flow satisfies the system ( |2.1D with the pressure determined by 

dk 1 



k and v are determined from (|3.9D . | 

Thus any pressure p may be obtained via a suitable potential r/. 

4 Hamiltonian Structure of Extremals 

We now briefly explore the Hamiltonian nature of (|2.1|) and the extremal 



equations (3^). If u is a smooth function of x and t, and h[u] is a function 



of X, t and the jet of u, let 

H[u] = I h[u]dx. 
Jr3 



Define 



5H[u] = / ( -^^,6u ) dx. 

We have the following result: 
Theorem 8 For the Hamiltonian (3.6) 



Thus the extremal equations (|3.2| ) may be written as 

Svr _ ^ 50 _ 5H 
'dt~~~66' 'dt~~6n' 



(4.1) 



These equations are canonical with respect to the natural symplectic form 
on L2(R3 : R^) X L2(R^ : R^) 



a;((Xi,yi),(X2,l2)) 



/ {Y2-X^-X2-Y^)dx 

yR3 



Thus we have expressed the extremal equations in terms of a canonical 
Hamiltonian system. 



5 Conclusions 



As described earlier, the Euler (impulse) equations ( |3.2| ) are not quite in the 
symmetric form that we obtain in the rigid body setting - i.e. we do not 
get a symmetric flow on two copies of the diffeomorphism group. However, 
it is possible to extend the analysis to this setting by factoring tt as 

= r o tp, (5-1) 

where 

^=voiP (5.2) 

and ijj evolves on the diffeomorphim group. Thus we do get symmetric 
equations for 4> and ^l' coupled to an interesting radial equation for r. This 
also has an analogue in the finite-dimensional setting - one allows P to be in 
Gl{n) and considers the polar decomposition P = RK where R is symmetric 
positive definite and K lies in SO(n). We shall describe the details of this 
analysis in a forthcoming publication. 

In addressing these issues, a deeper understanding of both the Hamil- 
tonian and variational structure as well as the geometry is needed. For 
example, we can obtain more insight into some of the calculations done in 
this paper as follows. Consider the Hodge projection P : X ^ X^oi tak- 
ing a vector field z to its divergence free part parallel to the boundary. As 
we have mentioned in the introduction, using the L2 pairing, this map is a 
Poisson map. Taking the L2 kinetic energy as the Hamiltonian on the un- 
constrained space X as well as on the constrained space Xvoi, we conclude 
that the corresponding Hamiltonian systems with their Lie Poisson bracket 
structures are mapped one to the other (including integral curves) by the 
Hodge projection. This simple remark is, in fact, the essense of what is going 
on in relaxing the divergence free constraints and in relating the Hamilto- 
nian structure in the formalism of Osledets, Buttke, and Kusmin. We have, 
in fact, shown some aspects of this remark in the above direct calculations. 
A deeper problem, to which we will return in other work, is to carry this out 
in material representation, where one needs a nonlinear Hodge decomposi- 
tion, similar to the Moser decomposition (a diffeomorphism group analogue 
of the polar decomposition) discussed in Ebin and Marsden [1970]. Many of 
these issues are addressed in work of Brenier; see, eg, Brenier [1999]. 

Acknowledgement: We would like to thank Peter Smereka for useful 
conversations. 
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